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Abstract 

Neutron stars are excellent emitters of gravitational waves. Squeezing matter beyond nuclear densities invites 
exotic physical processes, many of which violently transfer large amounts of mass at relativistic velocities, 
disrupting spacetime and generating copious quantities of gravitational radiation. I review mechanisms 
for generating gravitational waves with neutron stars. This includes gravitational waves from radio and 
millisecond pulsars, magnetars, accreting systems and newly born neutron stars, with mechanisms including 
magnetic and thermoelastic deformations, various stellar oscillation modes and core superfluid turbulence. I 
also focus on what physics can be learnt from a gravitational wave detection, and where additional research 
is required to fully understand the dominant physical processes at play. 
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1 Introduction 

The dawn of gravitational wave astronomy is one of 
the most anticipated scientific advances of the coming 
decade. The second generation, ground-based gravita¬ 
tional wave interferometers, Advanced LIGO (aLIGO 
Aasi et al. 2015b) and Virgo (Acernese et al. 2015), are 
due to start observing later in 2015 and in 2016, respec¬ 
tively. The first aLIGO observing run is expected to be 
a few times more sensitive than initial LIGO, with a 
full order of magnitude increase in strain sensitivity by 
^ 2019 (for a review, see Aasi et al. 2013). 

The inspiral and merger of compact binary systems 
(neutron stars and/or black holes) are commonly ex¬ 
pected to be the first detections with aLIGO. The most 
robust predictions for event rates come from obser¬ 
vations of the binary neutron star population within 
our galaxy, with an expected binary neutron-star de¬ 
tection rate for the aLIGO/Virgo network at full sen¬ 
sitivity between 0.4 and 400 per year (Abadie et al. 
2010b). But there are many other exciting astrophysi- 
cal and cosmological sources of gravitational waves (for 
a brief review, see Riles 2013). Loosely, these can be 
divided into four categories: compact binary coales¬ 
cences, bursts, stochastic backgrounds and continuous 
waves. Burst sources include gravitational waves gener¬ 
ated in nearby supernova explosions, magnetar flares, 
and cosmic string cusps. Stochastic backgrounds arise 
from the incoherent sum of sources throughout the Uni¬ 
verse, including from compact binary systems, rotat¬ 


ing neutron stars, and primordial perturbations dur¬ 
ing inflation. Gontinuous gravitational waves are almost 
monochromatic signals generated typically by rotating, 
non-axisymmetric neutron stars. 

The above laundry list of gravitational wave sources 
prominently features neutron stars in their many guises. 
While supranuclear densities, relativistic velocities and 
enormous magnetic fields are exactly what makes neu¬ 
tron stars amenable to emitting gravitational waves of 
sufficient amplitude to be detectable on Earth, it is also 
these qualities that makes it difficult to provide accu¬ 
rate predictions of the gravitational wave amplitudes, 
and hence detectability, of their signals. A positive grav¬ 
itational wave detection from a neutron star would en¬ 
gender great excitement, but it is the potential to un¬ 
derstand the interior structure of neutron stars that will 
make this field truly revolutionary. 

In this review, I provide a detailed overview of many 
proposed gravitational wave generation mechanisms in 
neutron stars, including state-of-the-art estimates of the 
gravitational wave detectability. These include gravita¬ 
tional waves generated from magnetic deformations in 
newly-born and older isolated radio pulsars (section 2), 
accreting systems (section 3), impulsive and continuous- 
wave emission from pulsar glitches (section 4), magne¬ 
tar flares (section 5), and from superfluid turbulence in 
the stellar cores (section 6). As well as discussing grav¬ 
itational wave detectability, I also concentrate on what 
physics can be learnt from a future, positive gravita¬ 
tional wave detection. 
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It is worth stressing that this is not a review of all 
gravitational wave sources in the audio band, but is 
instead designed to review the theory of gravitational 
wave sources from neutron stars. The article is but one 
in a series highlighting Australia’s contribution to grav¬ 
itational wave research (Howell 2015; Kerr 2015; Slag- 
molen 2015). 

2 Magnetic Deformations 

Spinning neutron stars possessing asymmetric deforma¬ 
tions emit gravitational waves. Such deformations can 
be generated through elastic strains in the crust (Bild- 
sten 1998; Ushomirsky et al. 2000), or strong magnetic 
fields in the core (Bonazzola & Gourgoulhon 1996). 
A triaxial body rotating about one of it’s principal 
moments of inertia will emit radiation at frequency 
/gw = where v is the star’s spin frequency. More 
precisely, a freely processing, axisymmetric body, with 
principal moment of inertia, /zz, and equatorial ellip- 
ticity, e, emits a characteristic gravitational wave strain 
(Zimmerman & Szedenits 1979) 

dTT^G /zz/gwE 

= 4.2 X 10-26 V—) \ 

V10“6/ \10ms/ \lkpc/ 

( 1 ) 

where d is the distance to the source. For comparison, 
the smallest, upper limit on stellar ellipticity for young 
neutron stars in supernova remnants comes from LIGO 
observations of Vela Jr. (G266.2-1.2), with e < 2.3 x 
10“^ (Aasi et al. 2014b). The overall ellipticity record- 
holder is for the millisecond pulsars J2124-3358 and 
J2129-5721 with e < 6.7 x lO"® and e < 6.8 x 10"®, re¬ 
spectively (Aasi et al. 2014a). Typical ellipticity con¬ 
straints for isolated radio pulsars are e < 10^ - 10® (Aasi 
et al. 2014a). 

2.1 Spin down limit 

An absolute upper limit on the gravitational wave strain 
from individual pulsars, known as the spin down limit, 
can be calculated assuming the observed loss of rota¬ 
tional energy is all going into gravitational radiation. 
The left hand panel of Figure 1 shows the current upper 
limits on the gravitational wave strain from a search for 
known pulsars (red stars; Aasi et al. 2014a), as well as 
the spin down limits for the known pulsars in the ATNF 
catalogue (Manchester et al. 2005)^. The observed grav¬ 
itational wave upper limits beat the spin down limit for 
both the Crab and Vela pulsars, while a further five pul¬ 
sars are within a factor of five. Also plotted in Figure 


^ http: / / WWW. atnf .csiro. au/people / pulsar / psrcat / 


1 are the projected strain sensitivities for aLIGO and 
ET (solid and dashed black curves respectively), and 
the strain sensitivity for the S5 run of initial LIGO (for 
details, see Aasi et al. 2014a). Age-based upper limits 
on the gravitational wave strain can also be calculated 
for neutron stars with unknown spin frequencies (Wette 
et al. 2008). 

2.2 Magnetic field-induced ellipticities 

As a general rule of thumb, a neutron star’s ellipticity 
scales with the square of the volume-averaged magnetic 
field inside the star, (B) (e.g., Haskell et al. 2008), im¬ 
plying the characteristic gravitational wave strain also 
scales as ho ^ (B)^. The dipole, poloidal component of 
the magnetic field at the surface of the star is inferred 
from the star’s spin period and it’s derivative, but very 
little is known about the interior field strength and/or 
configuration. 

A considerable body of work has been devoted to un¬ 
derstanding possible magnetic field configurations^, in 
part to answer the question of gravitational wave de¬ 
tectability, but also to understand the inverse problem 
regarding what can be learnt from future gravitational 
wave observations. Purely poloidal and purely toroidal 
fields are known to be unstable on dynamical timescales 
(e.g., Wright 1973; Braithwaite & Spruit 2006; Braith- 
waite 2007). Mixed fields, where a toroidal component 
threads the closed-field-line region of the poloidal held, 
commonly termed ‘twisted-torus’ helds, are generally 
believed to be dynamically stable (e.g., Braithwaite & 
Nordlund 2006; Ciolh et al. 2009; Akgiin et al. 2013), 
although it has been suggested this is dependent on the 
equation of state (Lander & Jones 2012; Mitchell et al. 
2015). Moreover, a recent study of the effect twisted- 
torus helds have on crust-core rotational coupling dur¬ 
ing neutron star spin down suggest these helds may 
be unstable on a spin down timescale (Glampedakis & 
Lasky 2015). 

While generating solutions to the equations of New¬ 
tonian and general relativistic magnetohydrostatics is 
a noble task, it is no substitute for understanding the 
initial-value problem that brings one towards realis¬ 
tic magnetic held conhgurations. In Newtonian simu¬ 
lations with soft equations of state (typically more ap¬ 
plicable to main-sequence stars than to neutron stars), 
Braithwaite and co. showed magnetic helds evolve either 
to axisymmetric, twisted-tori (Braithwaite & Nordlund 
2006; Braithwaite 2009), or non-axisymmetric conhg¬ 
urations (Braithwaite 2008), depending on the initial 
conditions. General relativistic magnetohydrodynamic 


^In Australia, work on possible magnetic field configurations in 
the cores of stars dates back to the 1960’s, where Monaghan 
(1965, 1966) calculated magnetic field equilibria in polytropes, 
which even included calculations of magnetic-field induced stel¬ 
lar ellipticities. 
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Figure 1. Left panel: Current upper limits on the gravitational wave strain from known pulsars (red stars; data from Aasi et al. 2014a) 
and the spin down limits for known pulsars in the ATNF catalogue (blue dots). Right panel: Gravitational wave strain predictions 
for known pulsars. The blue crosses and green squares are for normal neutron star matter with purely poloidal magnetic fields and 
A = 0.01, respectively [see equation (2) and surrounding text]. The red circles assume the neutron stars are colour-fiavour-locked phase 
[CFL; equation (3)] with {B) = lORp. In both figures, the solid and dashed black curves show the projected strain sensitivity for aLIGO 
and ET respectively, and the grey curve is the strain sensitivity for the initial S5 run assuming a two-year coherent integration (e.g., 
see Dupuis & Woan 2005) 


simulations ubiquitously show the development of non- 
axisymmetries (Kiuchi et al. 2011; Lasky et al. 2011, 
2012; Ciolfi et al. 2011; Ciolfi & Rezzolla 2012), however 
these simulations begin with a restricted set of initial 
conditions. The long-term evolution to stable equilibria 
in these systems is still an open question. 

The aforementioned uncertainty in possible mag¬ 
netic field configurations translates to uncertainty in 
gravitational wave predictions. The relative strengths 
of poloidal and toroidal components strongly affects 
the stellar ellipticity, and hence gravitational wave de¬ 
tectability (Haskell et al. 2008; Colaiuda et al. 2008; 
Ciolfi et al. 2009; Mastrano et al. 2011). Relatively stan¬ 
dard models suggest (e.g., Mastrano et al. 2011) 



where Bp is the poloidal component of the surface mag¬ 
netic field, A is the ratio of poloidal-to-total magnetic 
field energy (i.e., 0 < A < 1, with A = 0, 1 correspond¬ 
ing to purely toroidal and purely poloidal fields rep- 
sectively), and I have normalised to fiducial values of 
stellar mass and radius. Unfortunately, equation (2) is 
pessimistic for gravitational wave detection with aLIGO 
as the only neutron stars with Bp ~ 10^^ G are the mag- 
netars, but with spin periods of P ~ 1 - 10 s they emit 
gravitational waves at frequencies too low. The right 
hand panel of Figure 1 shows predictions for the gravi¬ 
tational wave strain for known pulsars from the ATNF 
pulsar catalogue (Manchester et al. 2005). The blue 
crosses assume a purely poloidal field i.e., equation (2) 
with A = 1, while the green squares assume a strong in¬ 


ternal toroidal component given by A = 0.01. The red 
circles are described below, but with normal neutron 
star matter the prospects for gravitational wave detec¬ 
tion from magnetic field-induced non-axisymmetries are 
grim. One therefore has to hope that Nature has been 
kind, and has provided neutron stars with strong inter¬ 
nal, toroidal components of the magnetic field. 


2.3 Generating large ellipticities 

Certainly, one expects newly bom neutron stars to have 
large internal toroidal fields and correspondingly large 
ellipticities; strong differential rotation combined with 
turbulent convection in the nascent neutron star drives 
an a-Q dynamo, which winds up a toroidal field as 
strong as ^ 10^® G (e.g., Duncan & Thompson 1992). 
Although the symmetry axis of the wound-up field is 
aligned to the spin axis of the star, the ellipsoidal star 
will evolve on a viscous dissipation timescale to become 
an orthogonal rotator through free-body precession, and 
hence optimal emitter of gravitational waves (Cutler 
2002; Stella et al. 2005; Dall’Osso et al. 2009). How 
such a strong field evolves over secular timescales is an 
open question. 

A recent spate of papers has shown that strong 
toroidal fields can also be achieved in systems in dy¬ 
namical equilibrium by prescribing different forms for 
the azimuthal currents in the star (Ciolfi & Rezzolla 
2013), and also by invoking stratified, two fluid stellar 
models (i.e., where the neutrons form a superfluid con¬ 
densate and the protons are either a normal fluid or 
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a superconductor; Glampedakis et al. 2012a; Lander 
2013). Interestingly, exotic states of matter in the core 
such as crystalline colour-superconductors allow for sig¬ 
nificantly higher ellipticities (Owen 2005; Haskell et al. 
2007). Recently, Glampedakis et al. (2012b) showed 
that, if the ground state of neutron star matter is a 
colour-superconductor, then the colour-magnetic vor¬ 
tex tension force leads to significantly larger mountains 
than for normal proton superconductors. They derived 
fiducial ellipticities for purely poloidal fields of^ 

PS 4.0 X 10-'(3) 

(4) 

where ^CFl respectively denote the elliptici¬ 

ties if matter is in a two-flavour phase (i.e., only the u 
and d quarks form superconducting pairs) and a colour- 
flavour-locked phase (where all three quark species are 
paired), and (B) is the volume averaged magnetic field 
in the core of the star. The red circles in Figure 1 
show gravitational wave predictions for a colour-flavour 
locked superconductor, equation (3), with the average 
internal field 10 times the observationally inferred sur¬ 
face field, i.e., (B) = lOHp. These results are fascinat¬ 
ing; a positive detection of gravitational waves from 
magnetic deformations in neutron stars is a fundamen¬ 
tal probe of the fundamental state of nuclear matter. 

It is worth mentioning that the discussion in this sec¬ 
tion pertains to the optimal case of a rigidly rotating, 
axisymmetric body rotating about one of it’s principal 
moments of inertia. Such a body emits monochromatic 
gravitational waves at a frequency of 2 i 4 . If the neu¬ 
tron star core contains a pinned superfluid, it will emit 
also at the spin frequency, v (Jones 2010). In general, 
a non-aligned rotator will emit at 2 z/, and a num¬ 
ber of frequencies straddling these values (Zimmerman 
1980; Jones & Andersson 2002; van den Broeck 2005; 
Lasky & Melatos 2013). For radio pulsars, such modu¬ 
lations would also be present in other observables such 
as pulse time-of-arrivals and radio polarisation (Jones 
& Andersson 2001; Jones & Andersson 2002). 

2.4 Oscillations in yonng neutron stars 

The strongest emitters of gravitational waves from mag¬ 
netic deformations are likely young, rapidly rotating 
neutron stars. Such stars may also emit through other 
channels, the most likely being the unstable r-mode, 
whose restoring force is the Goriolis force; see section 
3.2 for a detailed description of the physics of r-modes. 
Targeted gravitational wave searches of young neutron 


^It is worth noting that, in any kind of superconductor, the el- 
lipticity scales linearly with the volume-averaged magnetic field 
(Cutler 2002), cf. {B)'^ for normal matter. 


stars in supernova remnants can be adapted to set lim¬ 
its on the amplitude of such oscillations (Owen 2010). 
This was first done with a 12-day coherent search of 
S5 LIGO data targeting the neutron star in the super¬ 
nova remnant Gassiopeia A (Wette et al. 2008; Abadie 
et al. 2010a). The search has since been extended to nine 
young supernova remnants, with the most sensitive r- 
mode fractional amplitude being less than 4 x 10“® for 
Vela Jr. (Aasi et al. 2014b). Such a limit is encroach¬ 
ing ‘interesting values’ of the amplitude when compared 
to simulations that calculate the nonlinear saturation 
amplitude of various r-modes in young neutron stars 
(Bondarescu et al. 2009; Aasi et al. 2014b). 

3 Accreting systems 
3.1 Torque balance 

An observational conundrum drives research into gravi¬ 
tational wave emission from accreting systems; namely, 
the absence of accreting pulsars with spin frequencies 
14 > 700Hz (Chakrabarty et al. 2003; Patruno 2010). 
The argument is simple: measured accretion rates allow 
one to calculate the angular momentum being trans¬ 
ferred to the neutron star, which should be spun up 
to frequencies at, or near, the breakup frequency {v > 
1 kHz for most equations of state Gook et al. 1994). The 
common explanation is that the neutron stars are los¬ 
ing angular momentum through the emission of grav¬ 
itational radiation (Papaloizou & Pringle 1978; Wag¬ 
oner 1984; Bildsten 1998), although alternatives exist, 
most notably invoking interactions between the accre¬ 
tion disk and the companion star’s magnetic field (e.g.. 
White & Zhang 1997; Patruno et al. 2012). 

Balancing the accretion torques with gravitational 
wave emission allows one to estimate an upper limit on 
the gravitational wave strain independent of the emis¬ 
sion mechanism (Wagoner 1984): 



where i?io = i?/10km. Mi ,4 = M/1.4 Mq, Fx is the X- 
ray flux and F* = 10“® ergcm“^ s“^. The scaling 
arises only because of the spin-equilibrium assumption; 
slower rotators need larger ellipticities for torque bal¬ 
ance. If one builds a mountain without being in spin 
equilibrium, then the gravitational wave strain is sim¬ 
ply given by equation ( 1 ). 

Equation (5) implies the loudest gravitational wave 
emitters are the brightest X-ray sources such as the low- 
mass X-ray binaries, the brightest of which being Scor- 
pius X-1 (Sco X-1). Watts et al. (2008) utilised obser¬ 
vations of Sco X-1 and other known accreting systems 
to determine their gravitational wave detectability with 
current and future interferometers showing that, even 
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at torque balance, most systems will be very difficult to 
detect with aLIGO. 

The LIGO Scientific Collaboration has given periodic 
gravitational waves from Sco X-1 high priority, with 
multiple data analysis pipelines already fully developed 
(see Messenger et al. 2015, and references therein). The 
unknown spin period of the neutron star in Sco X -1 and 
the variable accretion rate that drives spin-wandering 
of the neutron star complicate these searches, but the 
best gravitational wave upper limit utilises the Side¬ 
band search (Sammut et al. 2014), which gives a 95% 
upper limit of ho ^8 x 10“^^ Hz at 150 Hz (Aasi et al. 
2015a). This is still above the torque balance limit for 
Sco X-1, « 3.5 X 10“^® (300Hz/^)^^^, but this is 

expected to be beaten with aLIGO observations (Whe¬ 
lan et al. 2015). 

Torque balance is an empirically derived limit that is 
independent of the mechanism generating the gravita¬ 
tional radiation. From a theoretical perspective, there 
are a number of ways in which sufficient energy can 
be lost to gravitational waves, including unstable oscil¬ 
lation modes (Andersson 1998), and non-axisymmetric 
deformations supported by the magnetic field (Cutler 
2002; Melatos & Payne 2005) or elastic crust (Bildsten 
1998; Ushomirsky et al. 2000). 

The torque balance limit for accreting systems with 
known spin periods is shown in the left hand panel of 
Figure 2 (blue dots), and for Sco X-1 (red curve) which 
has an unknown spin period. The LIGO, aLIGO and 
ET sensitivity curves in this figure assume two years 
integration; see Watts et al. (2008) for more detailed, 
realistic estimates of signal-to-noise ratios of these sys¬ 
tems. Note that current upper limits of gravitational 
wave emission from Sco X-1 utilise a semi-coherent, 10- 
day search, which therefore yields less stringent upper 
limits than the LIGO curve in Figure 2 (Aasi et al. 
2015a). 


3.2 Unstable oscillation modes 

Unstable oscillations modes, in particular r-modes, have 
received significant attention as potential sources of de¬ 
tectable gravitational waves. These inertial modes - 
toroidal modes of oscillation for which the restoring 
force is the Goriolis force - are generically unstable to 
a phenomena known as the Chandrasekhar-Friedman- 
Schutz (CFS) instability (Ghandrasekhar 1970; Fried¬ 
man & Schutz 1978). Consider an oscillation mode in 
a rotating star. The rotation induces both forward and 
backward propagating modes, but the backward prop¬ 
agating modes are being dragged forward by the star’s 
rotation. Although this retrograde mode continues to 
move backwards in the star’s rotating frame, if the star 
is rotating sufficiently fast it will be prograde in the 
inertial frame. The mode loses energy to gravitational 


waves, but carries with it positive angular momentum 
from the star. This positive angular momentum is sub¬ 
tracted from the negative angular momentum of the 
mode, which becomes more negative, growing the am¬ 
plitude of oscillation. As the mode grows, it emits more 
positive angular momentum, and continues to grow even 
faster. Hence, a sufficiently rapidly rotating star is un¬ 
stable to gravitational wave emission. 

At all stellar rotation rates, r-modes are retrograde 
in the comoving frame, but prograde in the inertial 
frame, implying they are generically unstable to the 
CFS instability (Andersson & Kokkotas 1998; Fried¬ 
man & Morsink 1998). A struggle therefore ensues be¬ 
tween the gravitational waves which drive up the os¬ 
cillation amplitude and the viscous damping that acts 
to suppress the mode amplitude. There exists a nar¬ 
row range in spin period and temperature for which 
gravitational wave emission dominates over viscous dis¬ 
sipation: at low temperatures, T < (few) x 10® K, vis¬ 
cous dissipation is dominated by shear viscosity, while 
at high temperatures, T > (few) x 10® K, bulk viscos¬ 
ity is the main culprit. Details of the r-mode instabil¬ 
ity window therefore depend sensitively on relatively 
unknown neutron star physics, including microphysics 
and complicated crust-core interactions (for a review, 
see Andersson & Kokkotas 2001). This somewhat sim¬ 
ple picture describing the r-mode instability window is 
inconsistent with observed spins and temperatures of 
low-mass X-ray binaries, implying the complete physi¬ 
cal picture behind this mechanism is ill-understood (Ho 
et al. 2011; Haskell et al. 2012). 

Strohmayer & Mahmoodifar (2014) recently analysed 
RXTE observations of the accreting millisecond pul¬ 
sar XTE J1751-305, finding evidence for a coherent 
oscillation mode during outburst that they attributed 
to either an r- or 5 -mode (for which buoyancy is the 
restoring force). The r-mode interpretation is the most 
interesting in the context of gravitational wave emis¬ 
sion; had the outburst occurred during the S5 LIGO 
run it would have been marginally observable (Ander¬ 
sson et al. 2014). However, doubt has been cast over 
the r-mode interpretation; Andersson et al. (2014) also 
showed that the mode amplitude required to interpret 
the observations necessarily leads to large spin-down of 
the neutron star, which is inconsistent with the neutron 
star’s observed spin evolution. 

3.3 Mountains 

On the other hand, permanent non-axisymmetries sup¬ 
ported magnetically or elastically may generate a sig¬ 
nificant gravitational wave signature. Accreted mat¬ 
ter accumulates on the neutron star, is buried, com¬ 
pressed, and undergoes a range of nuclear reactions 
(e.g., Haensel & Zdunik 1990). In the non-magnetic 
case, asymmetries in the accretion lead to compositional 


PASA (2015) 

doi: 10.1017/pas. 2015.XXX 


6 


Lasky 



^ 10-'' 

S' 

I 

N 9^ 

K 10"'^ 

^ jq-26 

'B 10-" 

CD 

^ TO 

" 10-'^ 

'i 10-'“^ 

-i-t 


gravitational wave frequency [Hz] 



-- 





111 

lllf ^ 


Figure 2. Left panel: Gravitational wave torque balance limit for known accreting millisecond pulsars and systems with burst 
oscillations assuming gravitational wave emission at twice the neutron star spin period (blue dots). Also shown is the torque balance 
limit for Sco X-1 (red curve) which has unknown spin period. Data is collated from Watts et al. (2008) and Haskell et al. (2015). Right 
panel: Gravitational wave predictions for magnetic mountains on known accreting X-ray pulsars, where the range is for magnetic field 
strengths at the onset of accretion of between = lO^'^ and 10^^ G (for details of the calculation, see Haskell et al. 2015). In both 
panes, the solid and dashed black curves show the projected strain sensitivity for aLIGO and ET respectively, and the grey curve is the 
strain sensitivity for the initial S5 run, assuming two years of coherent integration time. For comparison, current observational upper 
limits on Sco XI from LIGO are < 8 X 10“^® at 150 Hz, which utilise a 10-day, semi-coherent analysis (Aasi et al. 2015a). 


and heating asymmetries, which induce stellar defor¬ 
mations. Approximate expressions for the quadrupolar 
deformation can be derived in terms of the quadrupo¬ 
lar component of the temperature variation and reac¬ 
tion threshold energies (Ushomirsky et al. 2000). Such 
deformations are limited by the maximum stress the 
crust can sustain before breaking (Ushomirsky et al. 
2000; Haskell et al. 2006; Johnson-McDaniel & Owen 
2013), which typically gives larger gravitational wave 
estimates than torque balance (Haskell et al. 2015). 

Accretion also affects the structure of the neutron 
star’s magnetic field which, in turn, changes the ac¬ 
cretion dynamics. The magnetosphere funnels accreted 
matter onto the poles, at which point it spreads to¬ 
wards the equator, dragging, compressing, and burying 
the local stellar magnetic field (Hameury et al. 1983; 
Melatos & Phinney 2001). Such fields can lead to stel¬ 
lar deformations considerably larger than those due to 
the background magnetic field inferred from the exter¬ 
nal, dipole (Payne & Melatos 2004; Melatos & Payne 
2005; Vigelius & Melatos 2009), which may even gen¬ 
erate gravitational waves detectable by aLIGO or the 
Einstein Telescope (Priymak et al. 2011), particularly if 
the buried field is i? > 10^^ G (Haskell et al. 2015). This 
is shown in the right-hand panel of Figure 2, where the 
gravitational wave signal from accreting systems with 
known spin periods is shown assuming initial magnetic 
fields of i? = 10^° and 10^^ G (for details of this cal¬ 
culation, see Haskell et al. 2015). Finally, it is worth 
remarking that, given a positive gravitational wave de¬ 
tection, an in-principle measurement of cyclotron reso¬ 


nant scattering features in the X-ray spectrum may be 
able to discern magnetic or elastic mountains (Priymak 
et al. 2014; Haskell et al. 2015). 

4 Pulsar Glitches 

Pulsar glitches are sudden jumps in the neutron star 
spin frequency with wide-ranging fractional amplitudes, 
^ /S.v/v < 10“^ (Melatos et al. 2008; Espinoza 
et al. 2011; Espinoza et al. 2014). The exact mecha¬ 
nism driving a glitch is not fully understood, although 
it is clear that a sudden transfer of angular momentum 
occurs between the rapidly rotating superfluid interior 
and the outer crust (for a recent review, see Haskell & 
Melatos 2015). Such large angular momentum transfer 
lends itself to the generation of gravitational radiation 
through a variety of avenues, as both broadband burst 
emission from the glitch and as a continuous wave signal 
during the glitch recovery phase. 

4.1 Gravitational wave bursts 

Most theories of pulsar glitches rely on the general 
mechanism introduced by Anderson & Itoh (1975). The 
quantum mechanical nature of superfluids implies the 
neutron star’s rotation is attributed to an array of 
~ 10^® quantized superfluid vortices. These thread the 
entire star, but are pinned to lattice sites and/or crustal 
defects, and hence are restricted from moving outwards 
as the crust spins down through the usual electromag¬ 
netic torques. The superfluid core therefore retains a 
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higher angular velocity than the crust of the star, and 
a differential lag builds up between these two compo¬ 
nents. A glitch occurs when ^ 10^ - 10^® vortices catas¬ 
trophically unpin and move outwards, thereby rapidly 
transferring angular momentum to the crust. 

That so many vortices are required to unpin simulta¬ 
neously implies an avalanche trigger process must be at 
work. Such an avalanche is likely non-axisymmetric, and 
hence capable of emitting gravitational waves through a 
variety of channels. Warszawski & Melatos (2012) sim¬ 
ulated the motion of vortices, showing that a burst 
of gravitational waves emitted through the current 
quadrupole has a characteristic strain 


10 - ( 


Ar 


n 


10 - 


lO^rads ^ 

-1 / 7 \ -1 


1 kpc 


( 6 ) 


where Ar is the average distance travelled by a vor¬ 
tex during a glitch. The non-detection of gravitational 
waves from a glitch in the Vela pulsar during the fifth 
LIGO Science run in 2006 put an upper limit of /iq ^ 
10“^° (Abadie et al. 2011b), implying an upper limit of 
Ar 10“^ m and a lower limit on the glitch duration 
> 10“'^ ms. 

Figure 3 shows a histogram of estimates for the grav¬ 
itational wave strain using equation (6), where the 
glitches are taken from the ATNF glitch catalogue 
(Manchester et al. 2005)^ The vertical dashed line shows 
the empirical gravitational wave strain upper limit from 
the August 2006 glitch of the Vela pulsar (Abadie et al. 
2011 b). 


4.2 Glitch recovery 

The motion of vortices is expected to excite hydro- 
dynamic oscillation modes such as /- and p-modes, 
for which the pressure is the restoring force, and in¬ 
ertial r-modes (e.g., Andersson 1998; Andersson & 
Kokkotas 1998; Glampedakis & Andersson 2009). Al¬ 
though r-modes are expected to be the dominant emis¬ 
sion mechanism, recent estimates of their detectability 
with ground-based interferometers is pessimistic (Sidery 
et al. 2010). More optimistic estimates have been pro¬ 
vided for gravitational waves from meridional Ekman 
flows that are expected to couple the core and crust 
rotation post-glitch (van Eysden & Melatos 2008; Ben¬ 
nett et al. 2010), although it is worth noting the model 
is sensitive to many unknown, dimensionless parameters 
such as the Ekman number, normalised compressibility 
and Brunt-Vaisala frequency. Moreover, strong strati¬ 
fication may suppress Ekman flows to thin boundary 
layers in the outer regions of the star (Abney & Epstein 


^http: / / www.atnf.csiro.au/people / pulsar/psrcat / glitchTbl.html 
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Figure 3. Predicted gravitational wave amplitudes for known 
glitches using equation (6) with Ar = 10~^m (black, unfilled 
histogram) and Ar = 10“^ (red, filled histogram). The vertical, 
dashed line is the gravitational wave strain upper limit derived 
for the August 2006 glitch of the Vela pulsar during the LIGO S5 
run (Abadie et al. 2011b) 


1996; Melatos 2012), thereby reducing the gravitational 
wave strain estimates. Finally, it has recently been sug¬ 
gested that vortex avalanches that are believed to trig¬ 
ger glitches leave behind long-lived, large-scale inhomo¬ 
geneities on the vortex distribution that are potentially 
observable by aLIGO (Melatos et al. 2015). 

5 Magnetar Flares 

Perhaps the most exotic of all neutron stars are the 
magnetars; isolated neutron stars with external dipole 
magnetic fields Bp > 10^^ G. The decay of the immense 
magnetic field powers irregular bursts in hard X-rays 
and soft y-rays, with typical peak luminosities ~ 10^®— 
lO'^^ergs”^. Three giant flares have been observed in 
our Galaxy®, the largest of which was the extreme out¬ 
burst from SGR 1806-20 on 27 December 2004 which 
emitted a total, isotropic energy of 2 x 10"^® erg (Palmer 
et al. 2005). 

Giant flares are commonly attributed to catastrophic 
rearrangements of the magnetic field, either internal or 
external to the neutron star (e.g. Duncan & Thompson 
1992; Thompson & Duncan 1995). Strong coupling be¬ 
tween the magnetic field and the solid crust imply the 
latter will stress, and potentially rupture (although, see 
Levin & Lyutikov 2012), exciting oscillation modes in 
the star’s crust and core. Observations of quasi-periodic 


^Interestingly, LIGO non-detections of gravitational waves from 
nearby short-duration, hard-spectrum gamma-ray bursts GRB 
051103 (Abbott et al. 2008a) and GRB 070201 (Abadie et al. 
2012) rule out compact binary coalescences as their progeni¬ 
tors, implying these are most likely extragalactic magnetar giant 
flares. 
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oscillations (QPOs) in the tails of two giant flares (Is¬ 
rael et al. 2005; Strohmayer & Watts 2005) has engen¬ 
dered excitement about the potential new field of neu¬ 
tron star asteroseismology. In particular, the interpre¬ 
tation of QPOs as stellar magneto-elastic oscillations 
provide enticing potential to infer neutron star struc¬ 
tural parameters (see Levin 2006; Glampedakis et al. 
2006; Levin 2007; Sotani et al. 2007, 2008; Colaiuda 
et al. 2009; van Hoven & Levin 2011a,b; Levin & van 
Hoven 2011; Gabler et al. 2013, 2014; Huppenkothen 
et al. 2014c, and references therein). 

5.1 KiloHertz gravitational waves 

First attempts to understand gravitational-wave emis¬ 
sion from magnetar flares were optimistic. loka (2001) 
assumed an instantaneous change in the star’s moment 
of inertia from a rearrangement of the global, internal 
magnetic field can excite the neutron star’s fundamen¬ 
tal / mode (typically at frequencies, f ^ 1-2kHz), a 
study that was backed up by Gorsi & Owen (2011). 
These studies estimated that the energy in gravitational 
waves could be Sgw 10'^®-10"‘® erg, comparable to the 
energy emitted in electromagnetic waves, and certainly 
detectable in the Advanced Detector Era. 

In contrast. Levin & van Hoven (2011) showed that, 
while the order-of-magnitude estimates of loka (2001) 
and Gorsi & Owen (2011) are correct, in practice there is 
poor energy conversion into the /-mode, and prospects 
for gravitational wave detection in the Advanced Era 
are pessimistic. The analytic calculation of Levin & van 
Hoven (2011) looked at the energy conversion from a 
flare triggered in the star’s magnetosphere. Zink et al. 
(2012) and Giolfi & Rezzolla (2012) performed com¬ 
plementary, general relativistic magnetohydrodynamic 
simulations using their respective codes (Lasky et al. 
2011; Giolfi et al. 2011) for catastrophic reconfigurations 
of the internal magnetic field, also finding that the /- 
mode is not sufficiently excited to generate a detectable 
gravitational wave signature. Although these works had 
different predictions for the gravitational wave scaling 
as a function of magnetic field strength, they both pre¬ 
dicted approximately 10 orders of magnitude less en¬ 
ergy being emitted in gravitational waves than Gorsi & 
Owen (2011). 

Figure 4 shows a histogram of the possible gravita¬ 
tional wave energies emitted in the /-mode from known 
galactic magnetars® using the calculations of Giolfi & 
Rezzolla (2012) (shaded red histogram) and Zink et al. 
(2012); Lasky et al. (2012) (empty black histogram). In 
blue are the gravitational wave predictions from loka 
(2001) and Gorsi & Owen (2011), while the dashed black 
line gives the observed upper limit on the /-mode grav- 


®The known magnetars are taken from the McGill Magnetar cat¬ 
alogue (Olausen & Kaspi 2014). 



possible gravitational wave energy [erg] 


Figure 4. Predictions of the possible gravitational wave energy 
emitted in /-mode oscillations if giant flares were to go off in each 
observed galactic magnetar given the calculations of (Giolfi & 
Rezzolla 2012, shaded red histogram) and (Zink et al. 2012; Lasky 
et al. 2012, empty black histogram). Plotted in blue are the opti¬ 
mistic predictions of loka (2001) and Gorsi &; Owen (2011), and 
the dashed black line is the upper limit on the /-mode gravita¬ 
tional wave energy from the giant flare from SGR 1806-20 (Abadie 
et al. 2011a). 


itational wave energy from the 2004 giant flare in SGR 
1806-20 (Abadie et al. 2011a, see below). 


5.2 Low-frequency gravitational waves 

Although the /-mode couples well to gravitational ra¬ 
diation, one of the reasons it is not an optimal source 
for ground-based detectors is that it is damped through 
the emission of gravitational radiation in < 0.1s (De- 
tweiler 1975; McDermott et al. 1988). Longer lasting, 
lower frequency modes excited from magnetar flares 
have been suggested as potential sources of detectable 
gravitational waves (Kashiyama & loka 2011; Zink et al. 
2012; Lasky et al. 2012). To generate sufficient mass 
motions, and therefore a detectable gravitational wave 
signal, these modes need to be global core oscillations, 
and hence are unlikely the cause of the observed QPOs. 
A key uncertainty for these low frequency modes is their 
damping time, which underpins their gravitational wave 
emission. If these modes form an Alfven continuum of 
frequencies, where the QPOs arise as the edges or turn¬ 
ing points of the continuum (Levin 2007), then reso¬ 
nant absorption should quickly redistribute energy, re¬ 
sulting in short lifetimes (< 1 s) for the QPOs (Levin 
2007; Levin & van Hoven 2011; Gabler et al. 2011; van 
Hoven & Levin 2012; Huppenkothen et al. 2014c), and 
pessimistic estimates for gravitational wave emission. 
On the other hand, if global Alfven waves are excited 
in the core (Zink et al. 2012) then they may live for 
lOO’s of seconds, and be potentially detectable by third- 
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generation interferometers, or even aLIGO in the case 
of extremely strong (~ 10^® G) magnetic fields in the 
core (Glampedakis & Jones 2014). 


5.3 Current observational limits 

At the time of the giant flare from SGR 1806-20, the 
4 km Hanford detectors was the only interferometer op¬ 
erating. An upper limit of Agw < 7.7 x 10"^® erg was 
achieved around the 92.5 Hz QPO, with a significantly 
weaker constraint placed in the kHz range where one 
expects the /-mode. The most sensitive gravitational 
wave search of the giant flare was from the AURIGA 
bar detector (Baggio et al. 2005). They searched for 
exponentially decaying signals (with decay time 0.1s) 
in a small frequency band around 900 Hz. Subsequent 
LIGO/Virgo searches of regular flares from magnetars 
have yielded comparable gravitational wave limits to 
the original 2004 burst (Abbott et al. 2008b; Abadie 
et al. 2011a). 

While galactic giant flares only occur approximately 
once per decade, the regularity of magnetar flare storms 
implies they may provide an attractive alternative. If 
such flares also excite normal stellar modes, it is easy to 
imagine that stacking gravitational wave data from mul¬ 
tiple flares could allow for sufficient integration time. 
Indeed, such a gravitational wave search has been de¬ 
signed (Kalmus et al. 2009) and implemented to search 
for gravitational waves in LIGO/Virgo data (Abbott 
et al. 2009). These attempts have further been buoyed 
by the successful detection of QPOs in ordinary mag¬ 
netar flares (Huppenkothen et al. 2014a,b). 


6 Superfluid Turbulence 

Soon after the realisation that neutron star matter 
should form a superfluid condensate (Baym et al. 1969, 
and references therein), Greenstein (1970) suggested 
this superfluid should be in a highly turbulent state^. 
Such turbulence can be either hydrodynamic (e.g., Per¬ 
alta et al. 2005, 2006a,b, 2008) or quantum mechanical 
(e.g., Tsubota 2009; Andersson et al. 2003; Mastrano 
& Melatos 2005; Andersson et al. 2007; Link 2012a,b). 
Regardless, the turbulence is likely driven by differen¬ 
tial rotation between the superfluid core (with Reynolds 
number Re ^ 10^^) and the solid crust that is being 
slowly spun down by electromagnetic torques. 

6.1 Individual sources 

Although the turbulence is axisymmetric when av¬ 
eraged over long times, it is instantaneously non- 

^The beautifully succinct abstract is seldom seen today: ‘The 
neutron superfluid in most neutron stars should be in a highly 
turbulent state. If so, this turbulence drastically alters its rota¬ 
tional properties’ (Greenstein 1970). 
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axisymmetric, and therefore emits stochastic gravita¬ 
tional waves. The signal was first calculated in Per¬ 
alta et al. (2006a) by numerically solving the spheri¬ 
cal Couette-flow problem with Re < 10®. A subsequent 
analytic calculation (Melatos & Peralta 2010) made a 
more robust prediction of the gravitational wave am¬ 
plitude and spectrum, showing in particular that the 
peak of the gravitational wave signal occurs near the 
inverse of the turbulence decoherence timescale, Tc- For 
reasonable neutron star parameters, this decoherence 
timescale is 


26 


AH 


10 rad s 


ms. 


( 7 ) 


where AH is the difference in angular frequency between 
the crust and the core of the neutron star, which can 
potentially be related to the star’s angular frequency, 
H, through the Rossby number, Ro = AH/H. For many 
pulsars, this places the peak of the gravitational wave 
signal in LIGO’s most sensitive band, although it is 
worth noting there is considerable uncertainty in AH. 
Finally, the root-mean-square of the gravitational wave 
strain is 


h 


rms 


8 X 10"^® 


AH 


10 rads 


-1 



( 8 ) 


Figure 5 shows the peak gravitational wave predic¬ 
tions for the superfluid turbulence model for observed 
galactic pulsars and millisecond pulsars with Ro = 10“^ 
(black points) and 10“^ (red points). It is worth stress¬ 
ing that these values of the crust-core shear, AH, are 
most likely exaggerated; in the same paper Melatos & 
Peralta (2010) showed that Ro < 10“^ for typical mil¬ 
lisecond pulsars. Calculations of vortex (Seveso et al. 
2012) and flux-tube (Link 2003) pinning suggest AH 
should vary, but will generally be < 10“^. One there¬ 
fore concludes that turbulence excited in newly born 
protoneutron stars, or extremely close rapid rotators, 
are the only possible source for Advanced Era gravita¬ 
tional wave interferometers. The blue square and star in 
Figure 5 show the peak gravitational wave emission for 
hypothetical neutron stars situated at d = 10 pc with 
spin periods p = 3 ms and 10 ms, and Ro= 10“^ and 
10 “^, respectively, the former of which would be de¬ 
tectable by aLIGO. 


6.2 Stochastic background 

A cosmological population of neutron stars with tur¬ 
bulent cores produces a stochastic gravitational wave 
background that peaks in aLIGOs most sensitive fre¬ 
quency band (Lasky et al. 2013). Although there are 
large uncertainties in the expected amplitude of the 
background due to a lack of understanding of AH, 
the shape of the spectrum is relatively unique, with 
it being well-approximated by a piecewise power-law 
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Figure 5. Peak gravitational wave amplitudes from superfluid 
turbulence from galactic pulsars with Ro = AQ/Q, = 10“^ (black 
points) and 10“^ (red points). The blue star and square are hy¬ 
pothetical, nearby {d = 10 pc) rapid rotators with spin periods 
p = 3 ms and 10 ms, and Ro = 10“^ and 10~^, respectively. 


^gw(/gw) = f^a/gw> with a = 7 and a = -1 for /gw < 
fc and /gw > fc, respectively. Here, /c is the population- 
weighted average of I/tc, and Hgw(/gw) is the energy 
density in the gravitational wave background as a frac¬ 
tion of the closure energy density of the Universe. 

It is worth noting that turbulent convection in main 
sequence stars also produces a gravitational wave sig¬ 
nal (Bennett & Melatos 2014). Interestingly, the loudest 
gravitational wave signal detectable on Earth may come 
from the Sun as, for frequencies /gw ^ 3 x 10“^ Hz, the 
Earth lies in the Sun’s near zone, in which the wave 
strain scales significantly more steeply with distance, 
hrms oc d~^ (cf. (X d~^ in the far zone; Cutler & Lind- 
blom 1996; Polnarev et al. 2009). This gravitational 
wave signal is most relevant in the micro- to nHz regime, 
and therefore most applicable to pulsar timing experi¬ 
ments. However, the gravitational-wave scaling at low- 
frequencies is uncertain as Kolmogorov scaling for tur¬ 
bulence breaks down below /gw < 10“® Hz. 

7 Conclusions 

Ground-based gravitational wave interferometers are al¬ 
ready contributing to our understanding of interesting 
astrophysical phenomena. But non-detections will only 
progress the held so far; the hrst direct detection of 
gravitational radiation will herald a new scientihc held 
of study, and will allow new insight into the most exotic 
regions of our Universe. 

While the hrst detection with aLIGO is expected to 
be from the inspiral and merger of a compact binary 
system, there are many unknown, and ill-understood, 
mechanisms that can generate gravitational waves with 


signihcant amplitudes from isolated neutron stars. This 
review highlights many of those mechanisms where we 
have some understanding of the key, physical processes. 
It also highlights the vast uncertainty of many of these 
predictions, showing that they rely on knowledge of neu¬ 
tron star physics beyond current capabilities. But this is 
what makes the held so fascinating; a positive detection 
of gravitational waves from any of the mechanisms dis¬ 
cussed in this article will allow an unprecedented view 
into the heart of neutron stars, where some of the most 
exotic physics in the Universe takes place. 
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